Search for: All records

A class of two-fast, one-slow multiple timescale dynamical systems is considered that contains the system of ordinary differential equations obtained from seeking travelling-wave solutions to the FitzHugh-Nagumo equations in one space dimension. The question addressed is the mechanism by which a small-amplitude periodic orbit, created in a Hopf bifurcation, undergoes rapid amplitude growth in a small parameter interval, akin to a canard explosion. The presence of a saddle-focus structure around the slow manifold implies that a single periodic orbit undergoes a sequence of folds as the amplitude grows. An analysis is performed under some general hypotheses using a combination ideas from the theory of canard explosion and Shilnikov analysis. An asymptotic formula is obtained for the dependence of the parameter location of the folds on the singular parameter and parameters that control the saddle focus eigenvalues. The analysis is shown to agree with numerical results both for a synthetic normal-form example and the FitzHugh-Nagumo system. 

Abstract Systems of activator–inhibitor reaction–diffusion equations posed on an infinite line are studied using a variety of analytical and numerical methods. A canonical form is considered, which contains all known models with simple cubic autocatalytic nonlinearity and arbitrary constant and linear kinetics. Restricting attention to models that have a unique homogeneous equilibrium, this class includes the classical Schnakenberg and Brusselator models, as well as other systems proposed in the literature to model morphogenesis. Such models are known to feature Turing instability, when activator diffuses more slowly than inhibitor, leading to stable spatially periodic patterns. Conversely in the limit of small feed rates, semi-strong interaction asymptotic analysis shows existence of isolated spike-like patterns. This paper describes the broad bifurcation structures that connect these two regimes. A certain universal two-parameter state diagram is revealed in which the Turing bifurcation becomes sub-critical, leading to the onset of homoclinic snaking. This regime then morphs into the spike regime, with the outer-fold being predicted by the semi-strong asymptotics. A rescaling of parameters and field concentrations shows how this state diagram can be studied independently of the diffusion rates. Temporal dynamics is found to strongly depend on the diffusion ratio though. A Hopf bifurcation occurs along the branch of stable spikes, which is subcritical for small diffusion ratio, leading to collapse to the homogeneous state. As the diffusion ratio increases, this bifurcation typically becomes supercritical and interacts with the homoclinic snaking and also with a supercritical homogeneous Hopf bifurcation, leading to complex spatio-temporal dynamics. The details are worked out for a number of different models that fit the theory using a mixture of weakly nonlinear analysis, semi-strong asymptotics and different numerical continuation algorithms. 

Abstract The pedestrian-induced instability of the London Millennium Bridge is a widely used example of Kuramoto synchronisation. Yet, reviewing observational, experimental, and modelling evidence, we argue that increased coherence of pedestrians’ foot placement is a consequence of, not a cause of the instability. Instead, uncorrelated pedestrians produce positive feedback, through negative damping on average, that can initiate significant lateral bridge vibration over a wide range of natural frequencies. We present a simple general formula that quantifies this effect, and illustrate it through simulation of three mathematical models, including one with strong propensity for synchronisation. Despite subtle effects of gait strategies in determining precise instability thresholds, our results show that average negative damping is always the trigger. More broadly, we describe an alternative to Kuramoto theory for emergence of coherent oscillations in nature; collective contributions from incoherent agents need not cancel, but can provide positive feedback on average, leading to global limit-cycle motion.